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EXECUTIVE SUMMARY 


This thesis provides closed forms which descnbe the statistics of the angular 
distnbution in azimuth and elevation planes for geometrically based single bounce spheroid 
model. This spatial model can be generated by revolving an ellipse, which lies on the xz 
plane of a 3D Cartesian coordinate system, about the x-axis. This ellipse has semi-major 
axis a and semi-minor axis b, respectively. The base and the mobile stations are located in 
the xy plane with the mobile station being at the center of the spheroid and the base station 
at the ongin of the coordinate system. It is assumed that the mobile station is transmitting 
and the base station is receiving. The base and mobile stations are separated by a distance 
Dy 

The scattering spheroid model is applicable to a macrocell environment in which the 
received multipath signals at the base station are onginated after a scattering from the 
surrounding environment about the mobile station. Assuming that the scatterers are 
uniformly distributed with a constant scatterer density function in the volume V as defined 
by the geometry of the spheroid, we derive the marginal angle of arrival (AOA) probability 
density functions (p.d.f). in the azimuth and elevation plane. 

Using the derived closed forms for both p.d.f.s we generate plots to examine how 
the parameters a, b and D affect the angular distribution of the received waves. From the 


marginal AOA p.d.f. in the azimuth plane we conclude that: 


Ik. The higher the ratio @/D, the lower the probability for the received 
multipath signals at the base station to be confined to a small angular 
region centered about the axis which connects the base and the mobile 


stations. 


wy The angular distribution in the azimuth plane is independent of the length 
of the vertical dimension Bb of the spheroid. 
3 Comparing the results of the spheroid model for small values of b with 


those of the circular scattering model presented in [Ref. 3], it was 


observed that there exists a significant difference between them. 


Additionally, the effect of the a/D ratio on the r.m.s angle spread in azimuth plane is 
discussed. From the plot results of the r.m.s angle spread in azimuth plane versus a/D, 
we conclude that for values of a/D up to 0.5, there exists a linear relation. However, for 
higher values of a/D, the relationship changes to a nonlinear curve with a positive 


gradient. 


Regarding the marginal AOA p.d.f. in the elevation plane, the following 
conclusions were drawn from the plots: 
Jee Similar to the p.d.f in the azimuth plane, the higher the ratio b/D the lower 
the probability for the received multipath components at the base station to 
be restricted in a small angular spread around the axis, which connects the 
base and the mobile stations. 
i The r.m.s angle spread which also depends on the semi major axis a, has 


a linear dependence with respect to the b/D ratio. 
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I INTRODUCTION 


A. BACKGROUND 


The outdoor radio environment contains a large number of obstacles, such as 
buildings, cars, hills, etc, which obstruct, reflect and scatter the wireless radio signals. As 
a result, the transmitted signal arrives at the mobile station or at the base station from 
many different angles and times of arrival; this is known as multipath. Another common 
phenomenon in cellular systems is the co-channel interference, and it is created because 
in a certain geographic region(cluster) there are other subregions(cells) that utilize the 
same set of frequencies. Multipath fading and co-channel interference can increase 
significantly the bit error rate and they don’t allow high bit rate data services in wireless 
communications systems. The presence of multipath and the motion of the receiver or the 
transmitter causes the received signal envelope to be time-varying and experience a 
Doppler shift that depends upon the angle of .arrival of the incoming waves. Concepts 
such as angular spread, time delay spread, and Doppler spread have to be taken in 
consideration during a wireless communication system evaluation. In particular, they 
should be taken into consideration in the design of smart antennas used for diversity, 
beamforming and emitter localization applications. 

The multipath effect alone can be countered by applying antenna spatial diversity 
which implements an antenna array consisting of multi-elements. Spatial diversity is an 
arrangement whereby signals from multiple elements are combined to result in a signal 


that fades less rapidly than the individual element signals. Ideally this can be 


accomplished if received signals at the two elements are decorrelated. The degree of the 
decorrelation between the two received signals by the two elements can be expressed by 


the spatial cross-correlation coefficient which is a function of the joint p.d_f. fo.(9.8) 


describing the angular distribution of the incoming waves in the three-dimensional space 
at the receiver. 

The interference limited factor can not be overcome by using spatial diversity. 
Switched beam arrays can be used with some success if the angle of arrival of the 
interfering signal is not near the angle of arrival of the desired signal [Ref. 1]. Only the 
beamforming arrays have the potential to completely suppress the co-channel interference 
(if the number of elements exceeds the number of co-channel interferers) improving the 
wireless communication system performance. The performance of these type of smart 
antennas depends on the spatial characteristics of the wireless radio channel. The term 
“spatial characteristics of the radio channel” is used to mean the statistics of the received 
signal envelope in the spatial domain. These spatial characteristics may not always be 
available from measurements. In such a case one needs to devise geometric models to 
depict the scattering process in the spatial domain. Another area where knowledge of 
angular distribution is valuable is in the determination of the power spectral density [Ref. 
1]. Therefore the angular distribution of the received multipath components is important 
in evaluating the performance of radiowave link between the transmitter and the receiver 


in wireless communications. 


Spatial channel models have been created and applied in order to be able to derive 
the statistics of the received signal envelope for both macrocell and microcell 
environments. The term macrocell enviroment assumes that the scatterers surrounding 
the mobile station have the same height or are higher than the antenna of the mobile. 
This results in a received signal at the mobile station which arrives from all directions 
after bouncing from the surrounding scatterers. In contrast, the base station’s antenna is 
located higher than the surrounding scatterers. Therefore the received multipath signals 
at the base station onginate after scattering from the area near to the mobile station 
[Ref.2]. In microcell enviroment both base and mobile station are deployed in the same 
height as the surrounding scatterers. This implies that the scattering process takes place 
near to the base and mobile station. The characteristic difference between the types of 
environments is that the angular spread of the received signals at the base station 1s 
confined to a small region for macrocell environments. In the microcell enviroment this 
region tends to be larger. 

In Table 1, which is a modified version of a table from [Ref 2.], are listed some 


representative spatial channel models. 


SPATIAL CHANNEL | DESCRIPTION OF THE} APPLICATIONS 
MODEL MODEL 


1. Lee’s Model 






Uniform 








Zs Discrete 
Distribution 







3. Geometrically Based 
Circular Model 


4. Geommetrically Based 
Elliptical Model 
















5. Gaussian Wide Sense 
Stationery Uncorrellated 
Scattering (GWSSUS) 













Scatterers are evenly spaced 
on a circular ring about the 
mobile at a distance R. Each 


of the scatterrers represents 
the effect of many other 


scatterers and hence is 
referred to as an effective 
scatterer. The base station 1s 
deployed a distance D from 
the mobile station. 

A model similar to Lee’s 
model with the N scatterers 
evenly spaced over an 
AOA range. 





In this model it is assumed 
that the scatterers are 
uniformly deployed within a 
radius R about the mobile. 
The locations of the 
scatterers are defined by a 
Spatial scatterer density 
function. 

Scatterers are uniformly 
distributed in an_ ellipse 
where the base station and 
mobile station are the focii 
of the ellipse. 








N scatterers are grouped 
into clusters in space such 
that the delay differences 
within each cluster are not 
resolvable within the 
transmission signal BW. 


















Predicts correlation 
coefficient using a discrete 
angle of arnval(AOA) 
model. A modified version 
of this model can account 
for a Doppler _ shift, 
assuming that the scatterers 
are moving with an angular 
velocity. 








Predicts correlation 
coefficient using a discrete 
AOA — model. The 
correlation predicted by this 
model falls off more quickly 
than the correlation in Lee’s 
model. 
It is suitable for macrocell 
environments and predicts 
the joint time of 
arrival(TOA) and AOA and 
marginal density functions 
at both base station and 
mobile station. 








It can be determined the 
marginal AOA, TOA, the 





joint TOA and AOA density 


function. It is applied in 
microcell enviroment. 







Provides an _ analytical 
model for the array 
covariance matnx. 












6. Gaussian Angle of} Special case of  the| Provides an analytical 
Arrival (GAA) GWSSUS model with a|model for the array 
single cluster and angle of | covanance matrix. 


arrival statistics assumed to 
be Gaussian distributed 
about some nominal angle. 





Table 1. Description of some representative spatial channel models. “From Ref. [2].” 


Figure | through 4 demonstrate the first four spatial channel models, listed in 


Table 1. 


m Effective 
scatterers | 


Ly 


d-el Base station... 





Figure 1. Lee’s Model “From Ref. [2].” 
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Figure 3. Geommetrically Based Circular Model “From Ref.[2]”. 
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Figure 4. Geommetrically Based Elliptical Model From Ref.[2]”’. 


Every spatial channel listed in Table 1 has one common feature; the received 
signal at the base station is assumed to arrive from a two dimensional plane connecting 
the tips of the transmitter and the receiver. Hence the angle of arrival (AOA) includes 
only the azimuth angle information. 

In this thesis we consider a 3D spheroid model as shown in Figure 5 This 
spheroid can be generated by revolving an ellipse, which lies on xz plane, with semi- 
major axis and semi-minor axis a and J, respectively, about the x-axis. The base station 
and the mobile station are located in the xy plane with the mobile station at the center of 
the spheroid, and the base station at the origin. It is assumed that the mobile station is 
transmitting and the base station is receiving. The distance between the two terminals is 


D with the requirement that the semi-major axis of the spheroid be less than or equal to 


the distance D. In the model considered here it is obvious that the angle of arrival (AOA) 


of the received signals includes both the azimuth and elevation plane coordinates. 





Figure 5. The spheroid geometry. 


B. OBJECTIVE 


Assuming that the scatterers are uniformly distributed around the mobile station 
as defined by the geometry of the spheroid, the objective is to derive the marginal 
probability density function (p.d.f.) angles of arrival (AOA) in the azimuth and in the 
elevation planes at the base station. After extracting these quantities, we investigate how 
the angular spread in azimuth and elevation plane are related with the dimensions of the 
spheroid and the distance between the transmitter and the receiver, as well as how they 


are correlated with each other. 
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Il. PROBLEM FORMULATION 


A. INTRODUCTION 


In this chapter we consider a single bounce scattering spheroid model. The model 
can be used in a macrocell environment where the base station antenna height is 
relatively high. As a result there will be no scattering from scatterers near the base 
station. Figure 6, shows the spheroid scatterer density geometry and the notation used to 
derive the marginal p.d.f.s. As is shown in the Figure 6, the location of each scatterer is 


defined by the Cartesian coordinates (x, y,2Z) or the spherical coordinates (7,,¢,,0,), or 


(7, , Oy ; By. where £ represents the complementary angle of 9, in elevation plane. 





Figure 6. The spheroid scatterer density geometry. 


The method for the derivation of the joint TOA/AOA p.d.f. and marginal AOA p.d.f. for 
both azimuth and elevation planes is based on the approach presented in [Ref. 3]. In the 
Figure 6, it is assumed that scatterers are distributed randomly in the volume V of the 


spheroid according to the spatial scatterer density function f,, .(x%,y,z). It is further 


assumed that the waves arrive directly at the base station after undergoing scattering only 


once from within the volume V. Multiple scattering of waves is not considered in this 
thesis. 

Before proceeding with the determination of the p.d.f.s in the azimuth and the 
elevation plane we would like to present the formulation for the joint TOA/AOA p.d.f 
and marginal AOA p.d.f in the general case. In section B we start with the derivation of 
the joint TOA/AOA p.d.f.. In section C, we show formulation for the marginal AOA 
p.d.f.s. In the following sections D through G we derive the marginal AOA p.d.f in 


azimuth and elevation plane for the spheroidal model. 


B. DERIVATION OF THE JOINT TOA/AOA P.D.F.’S 


Figure 7, shows a plane wave impinging on a base station after scattering from an 


obstacle, represented by a point *S’ in space with the ordered rectangular triples (x, y, z). 


S4(%Y,Z) 


Tb 


BS 


*, 
*e 
2, 
e, 
'e 
2 
'e. 
% 
te 
te 





Figure 7. The 3D used coordinate system. 


The scatterer probability density function f, ,, - ge We z) can be also expressed 


terms of spherical coordinates systems (7%,@,,0, ) as follows. The equations relating the 


spherical coordinates to the Cartesian coordinates are: 
rh, =x? + y? +27 (1) 


0, = cos" 





o, =tan™'(~) (3) 


x =n, sin &, cos @, (4) 
= lp sin 0, sin D, (5) 
Le On 0, (6) 


where (X, y,Z) and (%,9,,0,) denote the location of the scatterer in Cartesian and 
spherical coordinates respectively. The joint p.d.f f cai (7,,0,.9, | can be found using 


the well known transformation [Ref. 4]: 


X= lp sin 0, COS @, 


Fry by.0, "> 9b») = I (1,94 .9b) Fe.y,2(% Ys ZY = 7 Sin B sin by (7) 
Z =n, COs, 


where J(7,,,,9,) is the determinant of the Jacobian transformation given by: 


Ox dx Ox 


dy oy a ae 
I(7,, 94,9, ) = 3 Bd 38, =r, sind, , 0<0, <7 (8) 


"Oca OcemeOr.. 
107, 9, AO, | 


Substituting the equation (8) into (7) we get for the scatterer probability density function 


oF 


b 


Ouch (7, QD, 0, ) with respect to spherical coordinates: 

hee ; 
Fir, 0, ("Ppp ) =m iN Oy fy y,2 (% Sin A, Cos M, ,m, sin, sin d, , 7, Cos, ). (9) 
If we want to express the scatterer probability density function dix ’, 0, (7,,0,.0, ) 


with respect to (7,@,,0,), where T denotes the time delay of the multipath component, 


we have first to find the relationship between 7% and the 7. Applying the law of cosines 


to the triangle BS,S,MS we get the following equation: 


15 


r? =r2 + D* —2r,Dcosd (10) 
where the cos 0 is called the ‘cosine direction’ given by: 
cosé = sin @, cos@,. (11) 


The delay corresponding to the total path traveled by the wave 1s: 


Gupte a ? . 
= btn 1, + afr? + D® + 2r,Dsin 8, cos, | (12) 


where ‘c’ is the speed of propagation of waves. 


Squaring both sides of the equation (12) and solving for 7 gives: 


1 ge 


b= 2(Dsin@, COs @, TC) ve) 


The joint TOA/AOA  p.df. ey P (7,0,,9, | can be found using the Jacobian 


transformation [Ref. 4] according to: 


Sn. dp On (1, Py Gp) 


= D*=-r"¢" 
Se Dp A, (7, by,O4) = |J (7% | 1s 2(Dsin 6, cosd,—7) (14) 





where J (r,,,. 0, } is the determinant of the Jacobian tranformation given by: 


wap 


__—-2(Dsin@, cos, -t¢)* (15) 


= : . 
D*c+t*c?—-21tc7 Dsin 0, COS, 


J I(%, Op) = 








Substituting the equation (15) into (14) we get for the joint TOA/AOA p.d-f. 
Ho, 6, G@ 0,9, ) 


D*c+t7c3-21c" Dsin @, cos@, 


Ft. 9096) = “T5 no, cose, 5% nates %-6) 16) 
Using the last equation (16) and the equations (9) and (13) we can express the joint 


TOA/AOA p.d.f. Heer yy (7,0,., } in terms of orginal scatterer probability density 


function ieee ye z] as follows: 


~17¢*)?(D*c+r°c3}-2%7Dsin 6; COs d ) 


- (7% S1n G, cos @, , %, Sin G, $in @ ,7%, cosG, )- 
8(D sin 6, cos ¢,-7)* fr.y.2% b Pp >% b SIN Dp 5 Fp ») 


Fc.0,.0, (T%pA) = sin O, 22 

(17) 

Note that the equation (17) expresses the joint TOA/AOA p.d.f. f oF) (z 0,0, ) as 

observed at the base station in terms of an arbitary scatterer probability density function 
fg, ow 

When the scatterers are uniformly distributed within an arbitrarily space A with a 


volume V, then the scatterer probability density function = ae (x, y, z) is given by: 


1 
7% ZE A 


fayelsona)| (18) 


Q, otherwise . 
If we consider this special case, the joint TOA/AOA p.d.f. f ae (z 0,9, ) can be 


written: 


(D*=17¢7)* (D*c+17c?-2%e7Dsin 8}, COS @, ) (19) 


T, 9p ,9,) = sin 
Sod, 0, ( Pp b) BAAD 8V(Dsin@, cos g,—tc)* 


where the range for the values of 7,0, ,0, has to be defined. 


C. DERIVATION OF MARGINAL AOA PROBABILITY DENSITY 
FUNCTION 


The marginal AOA p.d.f. with respect to spherical angle coordinates could be 


found by integrating the joint p.d_f. die One| (7,,0,,9, ) , aS given in equation (9) with 
respect to 7, over the range 7, (9,,8, ) Come (9,., ) , and then either with respect to 
0, (9, ) over the range O, (9, ) to 0, (9, ) or with respect @Q, (0, ) over the range 


D, (0, ) to —,, (0, ) depending on whether the marginal A.O.A p.d.f. is to be found in 


azimuth or elevation plane respectively. The terms 7 (%,4,) and r (@,,4,) are the 
lower and the upper limits of integration and depend on the azimuth angle @, and the 
elevation angle@,. The values of these limits can be found from the intersections of a 
straight line @, =constant, ¢~, =constant with the scatterer volume. The limits of 
integration in 6, , 6,,and@,.for the marginal AOA p.d.f. in azimuth plane can be 
defined as the smallest and the largest @,-values that bound a region R, which is the 
intersection between a plane at a specific @, value and the scatterer region volume. The 


above statement can be expressed in mathematical expressions as follows. 


The marginal AOA p.d_f. dios (9, ) in azimuth plane is given by: 


0 (9 y'% (9,.6, ) 
So, (6) = [og | fe0,.6, (r, 04.9, )dr,dO, > 
aa Ty) (9,.9, ) 


Py, (Y% ) to (9% ) 


- 


[ | sin, f,.y,.(% sin @, cos by, 7, sin O, sin g,, 7, cosO, )dr,dO, . (20) 
Ary (Op ) tin (94) 


The marginal AOA p.d_f. To, (0, ) in elevation plane is given by: 


(9,.8,) 
| Fie.0,.6, (7,0, 4, \dr, dd, = 
lp, (9,9, ) 


,, 
9,,(8,) 
d, (6, ) 


Op, (Op ) toy (Mp ) 

{|r sin, fx.y,-(% sin @, cos, ,7, sin, sin g,, 7, cos O, dr, dg, . (21) 
Oy (8; ) My (0, ., ) 
The limits for @, integration, which depend on the elevation angle 6, can be 


determined from the smallest and the largest azimuth angles that bound the region of 


intersection between the cone @, =constant and the scatterer volume. 


Equations (20) and (21) express the AOA marginal p.d-f.’s for any scatterer 


probability density function oe (exsyiz r When the scatterers are uniformly 
distributed in a volume with a constant probability density function She — (x, an the 


marginal AOA p.d.f. iil a) is given by: 


Py (Op) ty (Gp fo) 8» (% ) 
f,@)=t |  fr2sinO,dr,d0,= {lr2 @,.0,)- 73 (,.4, hinO,49, 
By, (8; ) ty (Py ) A, (%) 


(22) 
and the marginal AOA p.d.f. fo (6, ) is given by: 


Dy (O; )rp> (Hp Pb) oe » (9) 
fo, (@,) =+ | | r, sin6,dr,dd, = 4 fhe ( (0,9, )— 7, (0,9, sin 8, d¢,. 
Pry (5 ) ri (Gp ) = (9, ) 


(23) 


D. MARGINAL AOA P.D.F. IN AZIMUTH PLANE FOR THE SPHEROID 
SPATIAL CHANNEL MODEL 


Before we proceed to find the marginal AOA p.d.f. in the azimuth plane, we have 
to define the equation which descnbes the geometry of the spheroid in Cartesian and 


spherical coordinates. The spheroid in Cartesian coordinates is given by: 


(x-D)'+y? 22 
3 ++>=1 24 
"3 3 (24) 


and in spherical coordinates is given by: 


(a cos” 0, +b’ sin’ 6, )r,° — 2Db*r, sin, cos, + b*(D *(D? —a >) = Oo, (25) 


Instead of the polar angle 6,, its complementary 6 = ie: will be used hereafter for 
simplicity. Equation (25) can be expressed as: 

A(B)r, — B(B)r, cos ¢, +C =0 (26) 
where, A(f) =a" sin? B+b* cos? B, B(B)=2Db" cose = b?(D? _ a”) and 
A(B) = B°(B)cos’ ¢, —4A(B)C 

= 4b? |a? (a? sin* 8 + b* cos? B)- D’a’ sin? B — D*b’ cos” Bsin’ , | 
The two roots of the above equation are: 


a’ sin B+b7 cos” Bsin* Op 









Db? cos Bcos Pp - ab 






ry (%> B) = 


a Db? cos Bcosdptya a” sin’ B+b? cos” Bsin? Dp 


'b> (¢,, B)= 


(27) 
In terms of the angle f£, the marginal AOA p.d.f. So, (¢, ) in azimuth plane expressed 


as: 


rif 


fy, O) = 37 | ip (¢,.8)- 173 (0,. B)leos BiB (28) 


= fray (9, ) 





Lee ae ~ D* sin” @, 
a 


Be, SOL Te ao ee gr Rr it 
--silleess ha ealire -cicet lee cole, i, 5” pill desire Celta bea? 
~ewa4 3 


tr Sete hepeee ats entertain igmaaany h 


62, ei Sted 


o, =cons tan t 





> oS aCe 
Ja —D*° sin’ @, 


Figure 8. The intersected region R between a plane @, =constant 


and the spheroid. 


Figure 8 shows the limits of the angle P for a specific value of Q,, . The limits of 
integration in the elevation angle can be determined from the values of the angle / that 
bound the region R, which is produced by the intersection between the plane @, =constant 


and the spheroid. In our case the cross section of the spheroid with a plane @, =constant 


is an ellipse with semi-major axis equals to jae —~D* sin? g, and semi-minor axis 


equals to bila” — D? sin? @, and the center located at Dcos@,. For a general angle f, 


the straight line from the origin intersects the region R at two points 7, and 7, . At the 
two boundary points the roots 7, and 7, of the quadratic equation (26) in 7, coincide. 


This happens when the discriminant of the quadratic equation is set to zero. Carrying this 


out, the upper limit 4,,, is given by: 


(29) 





This equation is valid for @, inthe range of sin : (--2)< Q, <sin- (a). These 


correspond to the values of @, that bound the spheroid in the xy-plane. Taking all of 


the above into consideration, the marginal AOA p.d.f. in azimuth plane is given by: 


fame (9, ) 


fo, @)=4 Joos Bl (¢,.8)—73 (9,,B) HB, sin (-4)<¢,<sin"(¢). (30) 


The evaluation of the integral as defined in the equation (30) is straightforward. 


The following relations simplify some of the terms: 
ri (9,, 6) 7 r, (@,, B) - 


= (7, (0,.2) — ty (9, B)(7 (GB) i lb, (0,,P)N, (9,,8) 1 fs (Q,.2)). (31) 


r,, (0,,B)— 7, (9B) = (3 
r, (0, B) + %, (d,,B) = =A | (33) 

1s, y» B) + %;, (0,8) = ~*~ 26 | (34) 

r, (¢,,8)-%, (6, B) =a (35) 


Using equations (31) through (35) in equation (30) we get: 


aan (Opp 
2 2 2 " 
Ty, (Pp) == {(n, mle. TAG, 8 is Joos fdp = 
0 
Pow’) 153 Bycos? dy 4A BYC [are 4 om 
=3 RES pene a 8 A) COS Ops EC 
3V | ACB) ae ACh eos Pa . (36) 
0 


No 
to 


At a_ specific value of @ when f=£,,, the discriminant becomes 
A(B..,, )= B’(B.,,,) cos’ 9, —4A(B,,,, )C =0. Hence, writing 
A(B) = A(B)-A(B,.4. ) 
= [B’(B)cos? g, — 4A(B)C] - [B7( 6... Jcos” g — Alem eo): 
we See that 
A(B)= 4b? |a?(a? sin* 8 +b” cos” B)- D?a’ sin? B — D*b* cos” Bsin” ¢, | 
~ 4b? |a?(a? sin? Bra +b* cos” Bray IE D*a’ sin* Buax —D7b* cos” Borax Sin” o, | 
Cr), 
which gives: 
A=4b7|sin?_ - sin? Bfa*(D? — a?) +b7(a? - D’ sin’ ¢, ). (38) 


Using this the integral in equation (36) can be wnitten as follows: 





2 if aw 2 {fe 2 > Bmax 
p (¢ )= 16D7b* cos” ¢, 4 a7|D*-a* “—D* sin™ | Cos cost Bist Bre S008. Brox sin’ 6 oss SEE Eree eG 
d,\Tb/ av (a sin? B+b* cos” B)° 


_ 40-22? (Oa ED in)" Po eePsin Boa si0 ag, (39) 


3V (a*sin* B+b* cos” B)* 


The first integral in equation (39) can be simplified as: 


Tt ae 3 B|sin? Bay -sin? B “| " cos B(I-sin? ? B)afsin? Bmax —sin® B (40) 
Soiled elena GE, ae ee ep 


(a* sin? B+b" cos” B)? (a* sin? B+b7 cos” B)° 


On using the substitution t=sinf,and t,,, =sin Byox the last integral becomes: 


T cos B(1-sin* eis Bmax Siti” Pip a -TE (ee ee t,t? (41) 


(a” sin” B+b7 cos” B)° 


; ble PP 


Hence the term jean can be wnitten as follows: 


l-1° = <= ae +( +(a? -b?}?}) 
Using this in equation (41) we get: 

(1-1? io -{" mn max ~ Imax ~ 
eal pee os pee 


Furthermore, writing t=... sina, we get for the last two integrals: 


us 
max a 
| “= g VP na? =| I" max cos” 04 AQ (43) 
: (b?+(a?-b?)2) 4 (52+(a?—b?)12max sin? or) 


and 


tng : 

if t VP max? ee LA) ee, | = | t~ max cos’ @& 'o4 (44) 
= | — est iy 

A (p? +(a°—b 2} ‘ (b+(a2—b?)1? max sin? a) 


Similarly, the second integral in equation (39) can be expressed as 


x 
T cos By sin® Bax sin” B Bmax sin” B - fem? _ Vt manly Pale 17 4». COS” ee. (45) 
0 (p? Ge ane )sin~  (2+(a2-p?)sin2 By 0 (b? +(a? je 2} 0 (b? a —p? 7 max SiN 2 af 


The last integrals described in the equations (43) (44) and (45) are evaluated using the 


formulas 2.563-1, 2.563-2, and 2.563-3 from [Ref. 5]. 


The marginal p.d.f. in azimuth is then given by: 


Fo, (¢, )= in| »)- if (d, ). sin (--4)< ¢, <sin (2) 








a (46) 
21 n2_ 2 Ape Lee 4 ner L, .2.4 
w here I, (d, )= nD? COs > a D = 3a Ab = = = a a b (47) 
{7 (a0? 
and (48) 


2 2 


with Y =a’ sin” Bay +b* cos” Barry - 


Os MARGINAL AOA P.D.F IN AZIMUTH PLANE FOR THE SPECIAL 
CASE OF THE SEMI-MINOR AXIS OF THE SPHEROID GOING ZERO 


In this section we examine the special case of marginal AOA p.d-f., given by 
equation (46), for values of b—»0. In order to evaluate the p.d.f. for this case the 


following simplifications are made on equations (47) and (48): 
v=a’sin? B., +b” COS” Ba, = b7 cos” Bmax n + Stan? Borax} 
where tan G,,, iS given by equation (29). Since b-—>0 the last equation can be 


pep cos” . (49) 


expressed as follows: YS 22 
Ve! 


The terms vy ~b* and (y ~4?)(D? oe a? | can be expressed as: 


pi prog: 
=o =p and (y -b*)(D? ~a*)z b? (a? — D? sin? ¢,). (50) 


Plugging equations (47) and (48) to equation (46) the latter can be written as : 


a*\D?-a" |-b?|a*-D? sin? g, 


1\(6)- 1.6.) = 37°F) 
{- (p? - a? )y— 6?) + =F faa? ~ 4? i (4p4 ss 2a? y - a2b4 L| | 
(51) 
The term |(3a? — 46)? + (454 - 2026? y — a6! jt can be expressed as: 


(3a? — 4p? \y? - (4p4 ~2a?b? y - ab4 | =3a°y -2a*b* —a* —ab?y + 4b". (52) 


The first three terms in the nght side of equation (52) are second order terms, while the 
last two are fourth order terms with respect to b therefore can be neglected. Hence 


equation (52) reduces to: 





2 ee 2 
Bg Vata b qo ey 4b = | 20" on ip = _(p p?-a?]2+ D*=a" } 
= b= a’} D*~ cos” @ 
(53) 


Combining equations (49), (50) and (53) and taking also in consideration which terms are 


second or fourth order with respect to ‘b‘, equation (51) can be wnitten as follows: 


11(@, )— 12(@) = ele a” + D* sin” @, + D® cos” 0 Seep 2 po 


3V Deosg, a 7a D? cos, 
(54) 

Carrying out some simplifications in equation (54) we get: 
11 (4)- 1o(4) = 33 (a? - D? sin? 9, ). (55) 


Finally the marginal AOA p.d.f in azimuth plane of the spheroid for the values of b > 0 
is given by: 
So, (g, )= 32 (a? ~ D* sin 2 6.) sin (-.2)< ~, < sin!(2), (56) 


F. DERIVATION OF MARGINAL AOA PDF IN AZIMUTH PLANE 
THROUGH A SECOND APPROACH 


As we have mentioned before the marginal AOA p.d.f in azimuth plane is given 


by the: 
Bone Ths (9, .B) Bua. Ty (Pp » 
ie (O,)=F | [a sin Par, B =4 | ip COS Badr,dB . (37) 
Sores Tin (P% ?,.B) ‘ee in (9 B) 








b oe) 49 ’ 4 
— Ja*>-D* sin’ ¢, 
a 


Dcos¢, -Acos(z-@) 


Dcos 9, 


ja* —D* sin’ ¢, 


Figure 9. The 2D coordinate system used for the derivation of the azimuth p.d_-f. 


Figure 9, shows the intersected region between a plane at a specific @, value and 
the scatterer region volume. In the plane @,=constant, we adopt a 2D Cartesian 
coordinate system x;y;. The terms of the integral in the equation (57), y,cosf and 
r,dr,dP can be written as: 

r,cos 8B = Dcos@+Acosa@ and r,drdB = dx,dy, = dx,dy,. (58) 


In the x2y2 2D Cartesian coordinate system as shown in Figure 9, the term Acos@ = x,. 


According to the above equation (58) equation (57) can be written as: 


Bmax %2 YB) eya°—D* sin* % \a?-D* sin” 9, 
fy, Gp) = + ic cos Br, drdp = 4 | [(D cos, + Acos @ dx» dy, 
Sma "by (¢,,B) 0 0 


(39) 
The integral in equation (59) with respect to polar coordinate system A and @ can be 


expressed as: 


DA: 
f.()=+ | |(Doos¢, + Acosaydada. (60) 
a=0 0 


Taking in consideration that the term Acos@ is equal to x2, since x2 is an odd function, 


equation (60) can be written as: 


2 ax 
le Oar, F Tiweose, +X» WdAda=+ [| Pcosd, Addda (61) 
a=0 0 Ag, (» ) 
where A, (9, ) is the area of region R. 
Equation (61) is simplified as: 
f.,(G,) = + Deos¢, {[ Adada =F A, (¢,) (62) 


A(Q,) 
Where Ay, (9) =7 2 (a —D* sin* @, }. Substituting the above term in equation (62) we 
get for the marginal AOA p.d.f. in azimuth plane: 
Ig, (g,)= 3D e086 a? — D* sin ° O, sin (--2)< QO, S sin'(4). (63) 


The last solution approach shows that the azimuth angular distribution is independent of 
b and coincides with equation (56). Hence, we conclude that the limiting equation (56), 


valid for b—>0, is also true for all D. 


Finally, using the closed form expression for the p.d.f. we evaluate the r.m.s 


angular spread in closed form: 
Drs = 29 sr (64) 


where 


sin“!(2) 
tha = [5 f(s) = bin (g)f +2 [2 asin "Ug la—25)-284 3.05, 65 


“I2) 


nN 











to 


sin 


G. MARGINAL AOA P.D.F. INELEVATION PLANE FOR THE SPHEROID 
SPATIAL CHANNEL MODEL 

As noted in section C of the present chapter, the marginal AOA p.d.f. in elevation 
plane is given by: 


%, (8B) r, (9.8) 


fl B)=t [in cos fdd, = 37 ni ?,,B Vee (¢,, B)|cos Pdg, (66) 
?, (B) Ty (%.B) 9p, (B) 
Pp, (B)= max 
f(B)=% | & (9,.B)- 173 (9,.B)keos Bag, . |B < By (67) 
0 


The maximum elevation angle f,, is obtained from the values that bound the 
spheroid in the xz-plane and can be found analytically from the equation (29) by setting 


the azimuth angle @, equals to zero: 


(68) 





By = tan” (7s) ° 


Using the equations (31) through (35), the integral in equation (67) can be simplified as 


follows: 
Omax 
f(B)=s cos Br, = le Me +12 do, 
0 


ae 2058 i VB" (B)cos” %-4A(B)C B)cos” 2 9,-4A(B)C [22(a)eas*o )cos” %; Cc io 


ET,\¢ Dee Ge id) ee. (0) 


Vins 
= So cos” ?, B*(B)cos? ?, -4A(B)Cdé, a 
0 


nee 
ara j ve" )cos’ ¢, -4A(B)Cd@, . (69) 


Substituting the functional forms of A(#),B(f) and C from equation (26) we get: 


a Prax 
fo(B)= aera | 0" PV B*(B)cos? by -44(B)Cad,y ~ 
0 


26° |p" “=a * eos 6 *B2(B)cos? & — 4A(B)C 
3a? sin 2 B+b* cos” BY V Te B)cos Pp 4A B)Cd$, . (70) 


The upper limit for the azimuth angle @,,, , for a fixed £, is the largest value of angle 
g, that bound the region R , which is generated by the intersection of the cone 
f =constant with the spheroid. In general, the intersection will include all the points 
Ih S% S%,, wheren, and 7m, are given in equation (27). Once again, at the boundary 
points it yields y,=7,. Thus @,, is obtained by setting the discriminant 


A(@, ) = B°(B)cos* ¢, —4A(B)C_ equals to zero and solving with respect to @, : 





Ores (mre PC (ae) | 25) 
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In order to be able to solve the integral in equation (70) we have to simplify the 
expression for the discriminant A(@, ) = B*(B)cos’ g, -4A(6)C under the square root. 
This is done applying the same approach as in previous section D. The discriminant, at a 


specific value of { becomes equals to zero when the azimuth angle g, coincides with 
the maximum value of ¢, =@,,,. Thus writing the discriminant as 
A(O,) = A(9,)— A(Prx) 

=[B’(B)cos* 6, —4A(B)C]—-[B?(B)cos’ ¢,,,, —4A(B)C], becomes 
A(@, ) = 424 cos” Bcos” @, ~4p?(D? = a*|\a’ sin” B+b* cos” B) = 

- laD254 cos” Bcos” Pmax ~46?(D? = a®|a? sin? B+b7 cos” B} (72) 
which gives: 
A(@, )= A(d, )— A(@max ) = 4D7b* cos” B(sin* Omax — Sin? Op). (73) 


Substituting equation (73) into equation (70) gives: 


Grane 
16D°b° cos* B 2 we a 
fs Dagar cos? Gp (Sin? bnax — SiN” 9 )d0, - 


3(a7 sin? B+b* cos” B)*V 


Onax 
4/n,2 2 Z 
___4D6*(D?-a? )oos? 6 [| \(sin? Gnax sin? 9,40, (74) 
0 


The above two integrals may be rewritten as: 


Pmax Giiae 
[ cos?  y(sin? dnax —Sin? $,)dd, =Sin max | cos? ¢ ae sin? g,)dd, (75) 


0 0 


Ornax Pmax 
and [ J(sin? max —Sin? 4 )ddy =Sindmax f ,|(1- —al— sin” $5 )d6p (76) 
SIN Omax 
0 0 








These integrals can be evaluated in terms of complete elliptical integrals first and second 


a) 


kind. Indeed with the help of the formulas 2.595-1 and 2.592-2 of [Ref. 5], these 


integrals take the form: 


P max Soe. 3% 
; a) ” 1] uel mT ener ; _ SIN“ Omax—l : 
sing COS” @p, [} ae sin” @, eo. = —— E(Sin gmax )- 3 K (sing) 


(77) 





and 


P max 
sing f [lt-—2>—sin? 6, lig, = EGsin4g,,,,)- (sin?_, - 1)K (sin yy.) (78) 


where K(sin@g,,,) and E(sing,,.) denote the complete elliptical integral of the first and 





second kind respectively [Ref. 5]. 


Substituting the results from the equations (77) and (78) into equation (74) we get 


the final expression for the marginal AOA p.d.f., in elevation f B ( Bp ) in a closed form as 


follows: 


16D°*b° cos* 
fp (B) = Ey AT ( +sin? Ove (Sil @—== cos” Qe Sin) Oo ))|- 


_ —4064(D*-0? kos? 6 _[E (gin Pmax )~ COS Pmax K (sin Onax) BI = By (79) 


3(a* sin” B+b* cos” B)*V 


with By =can~{ = :} sin On = Sf{1- GSE)” and V =4na*b. 
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Il. NUMERICAL RESULTS AND DISCUSSION 


Having derived closed form expressions for the marginal AOA p.d.f. in azimuth 
and elevation planes, we now present in this chapter some numerical results for each 
marginal AOA p.d.f.. Recall that D is the distance between transmitter and receiver, a is 


the semi-major axis of the spheroid, and b is the semi-minor-axis of the spheroid. 


A. PLOTS FOR MARGINAL AOA PDF IN AZIMUTH PLANE 


In this section we consider the cases of “D’, ’a’ and ‘b’ as given in Table 2 and 


generate some plots: 





Table 2. Parameters of D, a and b used in plots for marginal AOA p.d.f. in the azimuth 
plane. 


Figure 10 has been plotted in order to check how the azimuth angular distributions 
of the spheroid and the circular scattering models are related to each other. The circular 
scattering model is presented in [Ref. 3]. The angular distribution for the circular model 


is given by equation (80): 
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D 


2Dcos@, J D? cos” ,-D* +a" - -I({_ a Sal F 
fo, (W) = na” sin "-$)s 0, ssin (5) (80) 


0, else 
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Figure 10. Plot of fg, (%, ) versus, for D=1000m, a=100m and b=0.001m. Also shows 
the AOA p.d.f. for a circular spatial model for D=1000m and a=100m. 


In Figure 10, we plot the p.d.f. for a/D=0.1 and b/D=10°. Comparing the plot of the 


spheroid model for small values of “b* with this of the circular model, we can observe that 


there exists a significant difference between them. 


In Figures 11 and 12 we plot the p.d.f. for a@/D=0.1, b=50m and for a/D=0.2, 


b=150m using the two different formulas for the angular distribution in the azimuth 


34 








| 
— pdf (eq(46) for 
Jr------ [ee ee ea a ee a/D=0.1 b=50m 
—— pdf eq(56) for 
I a/D=0.1 
6r------ Se ee ar ee ee ees es Le ay ae ee eee eee eo ees 
I | 
| | 
5+------ aS Ss Sy a eg pag ee 
I | 
“dD | | 
ea a eee (ae ee oe ee = She eee |r ae ee eee 
Se = 
— I { 
I 1 
3 So 6 SSS SQ eos Sea SB Oo te Ss oe fi as a : ea 
I t 
I ! 
Ata oT ! a Bie PS Mei oo 
t I 
| | 
Le = ==S>> © <= $6) =o Se | = - a 
! I 
| I 
0 
-15 -10 -5 0 5 10 iS 


Angle of Arrival in Admuth Plane(degrees),0,, 


Figure 11. Plot of fy, (¢, ) versus¢, for D=1000m, a=100m and b=50m 


(using the two different formulas). 


— pdf (eq(46) for 
a/D=0.2 b=150m 

—— pdf eq(56) for 
a/D=0.2 
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Figure 12. Plot of fy (9, ) versus ¢, for D=1000m, a=200m and b=150m 


(using the two different formulas). 
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As shown Figure 11 and 12, the two different formulas give identical results for the 
angular distribution in azimuth plane. Additionally, we conclude that the higher the ratio 
a/D, the lower the probability for the received signals at the base station to be restricted in 
a small angular region centered about the axis which connects the base and the mobile 
stations. 

Finally in Figure 13 we plot the r.m.s angle spread in azimuth plane as given in 


equation (64) versus the ratio a/D. 


60 
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Figure 13. Plot of the r.m.s angular ns; spread versus a/D. 


Figure 13 shows that for values up to 0.5, of the ratio a/D, the profile of the r.m.s angle 
spread is almost linear. Note that the r.m.s angle spread in azimuth plane is independent 


of b as proved with the derivation of the second formula for the p.d.f in azimuth plane. 
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B. PLOTS FOR MARGINAL AOA PDF IN ELEVATION PLANE 


Corresponding 


Figure 





Table 3. Parameters D, a, and 5b used in plots for marginal AOA p.d.f. in the 
elevation plane. 


Figure 14 is generated to demonstrate the angular distribution in elevation plane 
for small values of bD. 


oY 





Angle of Arrival in Elevation Plane(degrees),B 


Figure 14. Plot of fg(@) versus 8 for D=1000m, a=100m 


and b=Im.. 


As we can see form the Figure 14 the angular distribution in elevation plane for small 
values of 6 looks like a dirac function as it is expected. 

Figure 15 and 16 are plotted for demonstration purposes, giving information of 
how the changes of the ratio b/D for a given ‘a’ affect the angular distribution in elevation 


plane. 
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1000m, a=100m 


)versus 8 for D= 


gure 15. Plot of fp(B 


and b=30m.. 


Fi 


Figure 15 shows the the angular distribution in elevation plane for a ratio b/D=0.03 and 


a/D=0.1. 
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Figure 16. Plot of fg(B) versus 8 for D=1000m, a=100m 


and b=50m. 
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Figure 16 shows the p.d.f for a ratio b/D=0.05 and a/D=0.1 As seen, compared to Figures 
15 and 16, the higher the ratio b/D the lower the probability to receive multipath 


components from a small region {py at the base station. 


Figure 17 1s generated, as set with Figure 16, in order to show the effect of the 


length of the semi- major axis a on the angular distribution. 
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Figure 17. Plot of fg(@) versus @ for D=1000m, a=200m 
and b=50in meters. 


Figure 17 shows the p.d.f for a ratio b/D=0.05 and a/D=0.2. Comparing Figures16 and 
17, it is seen that p.d.f in elevation plane is almost independent of the length of semi- 
major axis of the spheroid for small! values of the ratio a/D. Figure 18, is generated to 


show the effect of higher values of the a/D ratio on the angular distribution. 
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Figure 18. Plot of fp) versus # for D=1000m, a=500m, 
and b=50m. 


Comparing Figure 17 and 18, it was observed that there exists a difference between the 
two angular distributions. 
Finally in Figure 19 we plot the r.m.s angular spread of the angle of arrival in 


elevation plane, as given by the equation (81), versus the ratio b/D for specific values of 


a/D. 


(81) 





The mean value of the angular distribution at the base station is zero, and the integral is 


evaluated numerically. 


4l 
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Figure 19. Plot of the r.m.s angular Bims spread versus b/D 
for various ratios a/D,for the spheroid scattering model. 
As shown in Figure 19, there exists a linear relation between the r.m.s angular spread and 
the ratio b/D independently of the ratio a/D. 
Eventually for validation purposes of the two derived, in this thesis, closed 
forms of the angular distributions, we plot in Figure 20 both p.d.f. with a almost equals to 


b (6 is not made exactly equal to a to avoid dealing with limiting process). 
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Figure 20. Plot of fy (¢,) versusd, and f B (6) versus £6 for D=1000m, 
a=100m and b=99.99m and for D=1000m, a=100m and b=100m respectively. 


Figure 20, shows the angular distribution in azimuth and elevation plane for 


=b/D=0.1. From the plot we can see identical results for both angular distributions, as 


expected due to the symmetry of the model. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


aXe CONCLUSIONS 


In this thesis we proposed a geometrically based single bounce scattering 3D 
model. It has been assumed that the scatterers are distributed about the mobile with a 
constant density function in a volume as defined by the geometry of the spheroid. The 
distance between the base station and the mobile station is D. The focus of the effort of 
this research was in the derivation of the mathematical expressions which described the 
spatial statistics properties of the received signals in the presence of the three-dimensional 
multipath scattering. 

After presenting a general approach for the derivation of the joint TOA/AOA and 
marginal AOA p.d.f.s in terms of an arbitary scatterer probability density function, we 
presented the derivation of the angular distribution in azimuth and elevation plane for the 
spheroidal model in closed expression forms. 

Making use of these quantities we studied the behavior of the angular spread for 
each p.d.f as function either of the spheroid dimensions or of the distance between the 
transmitter and the receiver. For either angular distribution the model results in a higher 
probability of arrival along the line of sight direction. Also, we observed that the higher 
the ratio of each spheroid’s dimension over the distance D the lower the probability that 
the multipath components are confined to a small region. An excellent agreement exists 
in the plots which are generated from the two different formulas of the p.d.f in azimuth 


plane. 
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Analytical and numerical calculations were performed for the r.m.s angular spread 
in azimuth and elevation planes respectively. The effect of the ratio a/D or b/D on the 


corresponding angular spread was investigated. 


B. RECOMMENDATIONS 


Further research is required to enhance this model. It will be a interesting to derive 
the marginal time of arrival p.d.f or to investigate the effect of the angular spread on the 
correlation observed between the received signals an antenna array. Another area of 
interest should be the validation of this model. This could be done comparing with 


measurements taken over similar macrocell environments. 
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